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Abstract 



A random-matrix theory is developed for the adiabatic response to an ex- 
ternal perturbation of the energy spectrum of a mesoscopic system. The 
basic assumption is that spectral correlations are governed by level repulsion. 
Following Dyson, the dependence of the energy levels on the perturbation 
parameter is modeled by a Brownian-motion process in a fictitious viscous 
fluid. A Fokker-Planck equation for the evolution of the distribution function 
is solved to yield the correlation of level densities at different energies and 
different parameter values. An approximate solution is obtained by asymp- 
totic expansion and an exact solution by mapping onto a free-fermion model. 
A generalization to multiple parameters is also considered, corresponding to 
Brownian motion in a fictitious world with multiple temporal dimensions. 
Complete agreement is obtained with microscopic theory. 

PACS numbers: 73.20.Dx, 05.40.-Fj, 05.45.-Fb, 71.25.-s 
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I. INTRODUCTION 



This is a theoretical investigation of the adiabatic response to an external perturbation of 
the energy spectrum of a complex quantum mechanical system. We consider a Hamiltonian 
'H{X) which depends on a parameter X. The X-dependence of the energy levels En{X) 
shown in Fig. 1 by way of example, is taken from a calculation of the hydrogen atom in 
a magnetic field.S Only levels with the same cylindrical symmetry are shown. A weak X- 
dependence of the mean density of states is removed by a rescaling of the energy. What 
remains is an irregular oscillation of as a function of X. Two levels which approach 
each other are repelled as X is increased further, leading to a sequence of avoided crossings 
at which the derivative En = dEn/dX changes sign. The average En is zero, averaged 
either over a range of X or over a range of n. The correlator of En{X) and Em{X') is non- 
zero for nearby levels n, m and for nearby parameters X, X', and serves as a quantitative 
characterization of how the system responds to an external perturbation. 

Our investigation was motivated by a remarkable universality of the parametric corre- 
lations discovered by Szafer and Altshuler.i They considered a disordered metallic particle 
with the topology of a ring, enclosing a magnetic flux (measured in units of h/e). The 
energy levels Ei{(f)) depend parametrically on (j). The dispersion is characterized by the 
"current density" 

j{E, 0) = S{E - E,{<P))j^E,{<P). (1.1) 

Szafer and Altshuler applied diagrammatic perturbation theory! to compute the correlation 
function 



C{8E, 64>) = j{E, 4>)j{E + + 50), (1.2) 

where the overline indicates an average over an ensemble of particles with different impurity 
configurations. The result was that the correlator C{dE, S(j)) becomes universal for 6E = 0, 



with P = 2 and X = Sep. Eq. ( |1.3|) is universal in the sense that it contains no microscopic 
parameters which characterize the particle, such as the diameter L, the mean level spacing 
A, the Fermi velocity v-p, or the mean free path /. It holds for {A/EcY^'^ <^ 50 <^ 1, where 
Ec — hvpl/L"^ is the Thouless energy. 

Eq. ( |1.3| ) was proven for the case that the randomness in the energy spectrum is due 
to scattering by randomly located impurities. Numerical simulations indicated that it ap- 
plies generically to chaotic systems, even if there is no disorder and all randomness comes 
from scattering at irregularly shaped boundaries.! (The average in that case is taken over 
E and 0.) Further work on disordered systems by Simons and Altshuler ,0 based on a non- 
perturbative "supersymmetry" formalism,! has shown that Eq. ( |1.3|) with P = 1 and X = 6U 
applies if the external perturbation is a spatially fluctuating electrostatic potential Us{r). 
(The function s(r) should vary smoothly on the scale of the electron wavelength, with van- 
ishing spatial average.) These analytical investigations assumed non-interacting electrons. 
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Recent numerical simulations of a Hubbard model by Faas, Simons, and Altshuleiil have 
shown that Eq. ( |1.3|) remains valid in the presence of electron-electron interactions. The 
correlator ( |1.3| ) thus provides a universal quantum mechanical characterization of the re- 
sponse of a chaotic system to an external magnetic or electric field. Such universality calls 
for a random-matrix theory of parametric correlations. It is the purpose of this paper to 
present such a theory. 

The basic principle of random-matrix theory (RMT) is that the spectral correlations 
are dominated by level repulsion.! Level repulsion is a direct consequence of the Jacobian 
ni<j \ Ei — Ej\^ associated with the transformation from the space of N x N Hermitian 
matrices 7i to the smaller space of N eigenvalues Ei. Level repulsion is universal in the 
sense that it is fully determined bv the symmetry class of the Hamiltonian ensemble. There 
exist just three symmetry classes,El characterized by the number /?= 1,2,4 of independent 
components of the matrix elements of 7i: /? = 1 in zero magnetic field (real 7i), /5 = 2 in 
non-zero field (complex Ti), and /3 = 4 for strong spin-orbit scattering in zero magnetic field 
(quaternion H). The three ensembles are called orthogonal (/3 = 1), unitary {(3 = 2), and 
symplectic (/3 = 4). 

The Wigner-Dyson theory of random matrices yields a level-density correlation function 
K{SE) which is universal for level separations 5E greater than the mean level spacing A.i 
The function K[6E) measures correlations between the level density 

N 

n{E,X) = Y.SiE-E,{X)) (1.4) 

i=l 

at different energies E and E + 6E, but at the same value of the external parameter X: 



K{6E) = n{E, X)n{E + 6E, X) - n{E, X)n{E + 6E, X). (1.5) 
The universal limiting form of K in the Wigner-Dyson theory is 



The universal correlator (|1.6|) was first obtained from RMT in the context of nuclear physics, 
and then applied to small metallic particles by Gorkov and Eliashberg.0 Much later, it was 
derived from a microscopic Hamiltonian by Efetovi and by Altshuler and Shklovskii.i The 
microscopic theory shows that Eq. (|1.6|) holds for a disordered metal in the energy range 
A <^ 6E -C Ec- Numerical simulations have established that the Wigner-Dyson theory 
applies generically to systems with chaotic classical orbits,0 and also that it remains valid 
in the presence of electron-electron interactions.lli 

The level-density correlation function (|1.6| ) is thus universal in the same sense as the 
parametric correlation function ( p..3[ ). This suggests that it should be possible to derive Eq. 
( |1.3| ) by some extension of the Wigner-Dyson theory to parameter-dependent Hamiltoni- 
ans TC{X). We will show that the Brownian-motion model used by Dysonlll to construct a 
parameter-dependent ensemble of random matrices, yields parametric correlations in agree- 
ment with the microscopic theory of Altshuler, Simons, and Szafer.i'@. 

The outline of this paper is as follows. In Sec. II we formulate the problem of a random- 
matrix theory of parametric correlations and define the mapping onto Dyson's Brownian- 
motion model. The correlation functions which we will calculate are summarized in Sec. HI. 
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In Sec. IV we present an asymptotic analysis which yields the correlation functions in the 
limit that the dimension of the Hamiltonian matrix goes to infinity. An exact result for 
the correlation functions in the Brownian-motion model for a special ensemble is given in 
Sec. V, and compared with the large- result of the previous section. In Sec. VI we extend 
the theory to parametric correlations involving multiple parameters. We conclude in Sec. 
VII, by comparing the results of random-matrix theory with the microscopic theory. 
The results of the asymptotic analysis were briefly announced in a recent Letter.0 



II. BROWNIAN-MOTION MODEL 

Starting point of our analysis is Dyson's Brownian-motion model0 for the evolution of 
an ensemble of N x N Hermitian matrices as a function of an external parameter r. Dyson's 
idea was to regard r as a fictitious "time" , and to model the r-dependence of the distribution 
of eigenvalues P({i?„},r) by the one-dimensional Brownian motion of A^ classical particles 
at positions -E'i(t), £'2(1"), .. . En{t), in a fictitious viscous fluid with friction coefficient 7 and 
temperature (3^^. Level repulsion is accounted for by the interaction potential —ln\E — E'\ 
between particles at E and E'. The particles move in a confining potential V{E), which is 
determined by the density of states. 

With these definitions, P({-E'„},r) evolves according to the Fokker-Planck equationEl 

dP d ( dW ,dP\ , , 

»'({£„}) = - J^lnlEi - E,\+Y^V{Ei). (2.2) 

Eq. has the r — >■ cxd ( "equilibrium" ) solution 

Peq({K}) = ^-^e-^^, (2.3) 

where Z is such that Pgq is normalized to unity. Eq. ( p.3|) , for /3 = 1, 2, and 4, is the 
eigenvalue distribution in the orthogonal, unitary, and symplectic ensemble. i It has the form 
of a Gibbs distribution, with the symmetry index j3 playing the role of inverse temperature. 
The fictitious energy W contains a logarithmic repulsive interaction plus a confining potential 
V . The function V{E) is chosen such that Peq yields the required average eigenvalue density 
(which depends on microscopic parameters, but is assumed to be independent of r). The 
logarithmic interaction has a fundamental geometric origin: The factor exp(/5X]i<j In l-Ej — 
= \Ei — Ej\^ is the Jacobian associated with the transformation from the space of 
Hermitian matrices Ti. to the smaller space of eigenvalues En- 

The A^-dimensional Fokker-Planck equation ( p.l|) is equivalent to A^ coupled Langevin 
equations. 



The random force is a Gaussian white noise of zero mean, J-'i^r) = 0, and variance 
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r,{T)r,{T') = ^6,,6{r-T'). (2.5) 



The Fokker-Planck equation ( p.l|) and the Langevin equations (p.4|) are equivalent levels of 
description of the Brownian motion.lli 

The fictitious time r needs still to be related to the perturbation parameter X in the 
Hamiltonian 'H{X) of the physical system one is modeling. Furthermore, we need a micro- 



scopic interpretation of the coefficient 7. These issues were not addressed in Ref. |I3|, but 
are crucial for our purpose. Let r = coincide with X = 0, so that 



N 



P{{En}.Q) = \[5{E,-E% (2.6) 



1=1 



with the eigenvalues of 7i(0). For r > we then identify 

r = (2.7) 

This is the simplest relation between r and X which is consistent with the average initial 
rate of change of the energy levels: On the one hand, 



(E,(X) - E^f = X' ' + 0{X') (2.8) 

is of order X^ for small X, while on the other hand the ensemble average 

2r 



(Hr) - Efy = — + 0{t') (2.9) 

is of order r for small r, according to Eqs. ( p.4|) and ( p.5|) . The identification (|2.7|) also 
implies the relation 



2 fdE, 



\dX 



(2.10) 



between the friction coefficient and the mean-square rate of change of the energy levels. 

Eq. ( |2.1| ) or ( ^.41 ) is the simplest description of the Brownian motion of the energy levels 
which is consistent with the equilibrium distribution ( p.3|) . It is not the most general de- 
scription: 1. One could include the "velocities" dEn/dr as independent stochastic variables, 
and work with a 2X-dimensional evolution equation. In the case of Brownian motion in a 
physical fluid, the appropriate evolution equation is Kramer's equation.lli It describes the 
dynamics of a Brownian particle on the time scale of the collisions with the fluid molecules. 
Since the viscous fluid in Dyson's Brownian-motion model is fictitious, it is not clear what 
the appropriate 2X-dimensional evolution equation should be in this case. 2. One could let 
7 be a matrix function 7jj ({£"„}) of the configuration of energy levels. Such a configuration 
dependence (known in fiuids as hydrodynamic interaction) would be an additional source of 
correlations, which is ignored. That is the basic assumption of Dyson's Brownian-motion 
model, that the spectral correlations are dominated by the fundamental geometric effect of 
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level repulsion. The Brownian-motion model is known to provide a rigorous description of 
the transition between random-matrix ensembles of different symmetry.^ However, there 
exists no derivation of Eq. ( p.l| ) or ( |2.4| ) from a microscopic Hamiltonian. Here we apply 
the Brownian-motion model to fluctuations around equilibrium in the random-matrix en- 
sembles ( |2.3| ), and show that there is a complete agreement with the microscopic theory for 
disordered metals.i'§ 



III. CORRELATION FUNCTIONS 

We consider observables A{X) of the form 

N 

AiX)=J2aiE,iX)). (3.1) 

i=l 

A quantity of the form ( |3.1| ) is called a linear statistic on the eigenvalues of 7i(X). The word 
"linear" indicates that A does not contain products of different eigenvalues, but the function 
a{E) may well depend non-linearly on E. We assume that a varies smoothly on the scale of 
the mean level spacing A. [In particular, this excludes the case of a step function a{E).] The 
correlator of A at two parameter values X and X' is 6A{X)6A{X'), where SA = A — A. The 
overline denotes an average over a range of X at constant 6X = X' — X (or, alternatively, 
over an ensemble of statistically equivalent systems). Of particular interest is the integrated 
correlator 

/•oo 

Xa= d5X5A{X)5A{X + 5X). (3.2) 

JO 

To compute the correlator of an arbitrary linear statistic we need the density correlation 
function 



(3.3) 



S{E, X, E\ X') =Y.KE- E,{XME'- E^{X')) 

- (r^w^^E^mj [t.s{e' - E,{x')) 

The correlator of A at X and X' then follows from a double integration, 

/oo roo 
dE / dE'a{E)a{E')S{E,X,E',X'). (3.4) 
-oo J—oo 



We will also consider the correlator A{X)A{X') of the derivative A = dA/dX of the 



linear statistic ( |3.1| ). This correlator follows from the density correlation function 5* by 

/oo roo Pp' 

JE J_JE'a{E)a{E')^j^S{E,X,E',X'). (3.5) 

Alternatively, we can compute the correlator of A from the current correlation function 



6 



CiE,X,E',X') = J2E,{X)Ej{X')5{E - E,{X))5{E'- Ej{X')). (3.6) 



Since 



^ • d 

A{X)=Y.E.{X)—a{E,{X)l (3.7) 

•i=i * 



one has, upon partial integration, 



— ; TOO roo fi'^ 

A{X)A{X') = J JEj JE' a{E)a{E')^^C{E, X, E', X'). (3.8) 



Comparison of Eqs. ( |3.5|) and ( |3.8|) shows that the density and current correlation functions 
S and C are related by 

g§^C(E, X, E', X') = g^S(E, X, e, X'). (3.9) 

We assume that S{E, X, E', X + 6X) = S{E, E', 6X) depends only on the parameter in- 
crement SX. When considering a particular physical system, such as the hydrogen atom in 
a magnetic field, this may require a rescaling of the energy levels, to eliminate a systematic 
drift in Ei versus X (cf. Fig. 1). S ince, by definition , S{E, E', 5X) = S{E', E, -5X), the cor- 
relators 6A{X)6A{X + 6X) and A{X)A{X + 6X) are even functions of 6X. Furthermore, 
we have the sum rule 

roo — : 

/ dSX A{X)A{X + 6X) = 0, (3.10) 
Jo 

for any linear statistic. 

In the following sections we will compute the density correlation function S{E,E',X) 
from the Brownian-motion model described in Sec. 0. In view of Eq. (|3.3|) and the identifi- 
cation (p.7|), we have the relation 



/oo roo roo roo I \ 

dEl ■■■ dE% dEi--- dEN\Y.^{E- E^)5{E' - Ej) 
~oo J —oo J —oo J —oo \ij I 

X P.^{{El]) (P({i?J,X2) -Pe,({i?n})) (3.11) 

between the density correlation function and the solution P({i?„},r) of the Fokker-Planck 
equation ( p^.l| ) with initial condition ( |2.(j| ). Once we have S", the current correlation function 
C and the correlators of A and A follow from Eqs. ( |3.4| ), (p.5|), and (p.9|). 



IV. ASYMPTOTIC SOLUTION 

In this section we compute the large-iV asymptotic limits of the density and current 
correlation functions S{E, E', X) and C{E, E' , X). By "asymptotic" we mean that the 
expressions obtained hold in the limit ^ cxd in the energy range \E — E'\ ^ A for all 
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X and in the parameter range X ^ for all E' . A justification of our asymptotic 

analysis will be given in Sec. 0, when we compare with an exact result for (3 = 2. We assume 
that the N ^ 00 limit is accompanied by a rescaling of the confining potential V{E) in Eq. 
( p. 21) , such that the mean density of states remains the same. An explicit example of such 
a rescaling is given in Sec. |V[ 

The first step in the analysis is to reduce the Fokker-Planck equation ( p.l|) to an evolution 
equation for the average density of eigenvalues 

/oo roo ^ 

dE^--- dENP{{En},r)Y,5{E~E{). (4.1) 
-00 J —00 

This problem was solved by Dysonlll in the limit A^ — > 00, with the result 

d 



p{E, t) — i V{E) - J JE' p{E', r)ln\E-E'\ 



(4.2) 



Corrections to Eq. ( ^4.2| ) are smaller by an order A^ ^ In A^. To the same order, the equilibrium 
density Peq{E) (defined as in Eq. ( |4.1| ) with P replaced by Pcq) satisfiesEl 

^i^ViE)- J°^JE'p,^{E')\n\E-E'\^=0. (4.3) 

To make this paper selfcontained, we present Dyson's derivation of Eq. ( [4. 21) in the Appendix. 

The next step is to reduce Eq. ( |4.2| ) to a diffusion equation by linearizing p around peq- 
This is consistent with the large-A^ limit, since p is of order A^ while fluctuations in the 
density are of order one.0 We write p{E, r) = PeqiE) + 5p{E, r) and find, to first order in 
6p, 

D{E, E') = -7- 111 |B - (4.5) 

Eq. ( [4.4| ) has the form of a non-local diffusion equation, with diffusion kernel ( [4.5| ). 

To proceed we assume a constant density of states Peq{E) = po = 1/A over the en- 
ergy range of interest (which is the energy range where the function a{E) in the linear 
statistic (3J-) differs appreciably from zero). The diffusion kernel can then be taken to be 
translationally invariant, D{E,E') = D{E' — E), with Fourier transform 

D{k)= r dEe"'^D{E) = (4.6) 
j-00 7|/c| 

Equation ([4.4|) becomes an ordinary differential equation in fc-space, with solution 

Sp{k, r) = 6p{k, 0) exp[-k'^D{k)T]. (4.7) 
In view of Eq. (|2.6| ) , the initial condition on the eigenvalue density is 

N 

p{E,0) = Y.S{E-E^). (4.8) 

i=l 
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We define the equilibrium average (/)eq of an arbitrary function of the initial 

configuration by 

/CX) fOO 
dEl-. dE%P^^{{El})f{{El}). (4.9) 
-oo J — oo 

Using also definition ([4.1|), Eq. ( ^.111 ) for the density correlation function S{E, E' , X) can 
be written as 

S{E, E', X) = {p{E, 0)p{E', X2))eq - Peq(i?)Peq(i?') 

= {Sp{E,0)Sp{E',X')),^. (4.10) 

In the second equality we have used that (p(£',r))eq = Pcq{E). The correlation function 
K{E, E') is defined by [cf. Eq. (|T3|)] 

K{E, E') = -{6p{E, 0)6p{E\ 0))eq = -S{E, E', 0). (4.11) 

Over the energy range of a constant density of states, the correlation functions S{E, E', X) = 
S{E' — E,X) and K{E,E') = K{E' — E) are translationally invariant, with Fourier trans- 
forms S{k,X) and K{k). According to Eqs. (^, ( ^TD] ), and ( PH ), we have 

S{k, X) = -K{k) exp[-k'^D{k)X'^]. (4.12) 

The function K{k) is known.i In the limit oo, one has asymptotically 

K{k) = (4.13) 
up 

independent of V{E)B Eq. ( ^ ) is the Fourier transform of Eq. ( |1.6| ) , and holds for energy 
scales k^^ ^ A large compared to the mean level spacing. (This is the relevant regime, 
since the function a{E) in the linear statistic ( p.l|) is assumed to be smooth on the scale of 
the level spacing.) 

Combining Eqs. ( [4. 61) , ( [4.12|) , and (|4.13| ), we conclude that the density correlation func- 
tion is given by 

S{k,X) = ^eM-e\k\), (4.14) 
up 

i^X{^p,hf'\ (4.15) 
The E'-space correlation function becomes, upon inverse Fourier transformation. 

The current correlation function C{E, E' , X, X') = C{E' — E,X' — X) is obtained from S 
by means of relation ( p.9| ), which in /c-space takes the form 

C{k,X) = -^^^S{k,X). (4.17) 
9 



We find from Eqs. ( [4. 141) and (|4.17|) : 



C{k,X) = '^-^{l-2e\k\)eM-e\k\). 

PI 



C{E,X) 



ln(e^ + 



(4.18) 
(4.19) 



The asymptotic results for the correlation functions given above can be used to compute 
the N ^ oo limit of the integrated correlator x^, defined by Eqs. (|3.2| ) and (O). The 
fc-space expression for Xa 



1 

2nJo 



dX 



dk\a{k)\^S{k,X), 



a{k) 



dEe'''^a{E). 



(4.20) 
(4.21) 



Substituting the asymptotic formula ([4.14|) , and carrying out the integral over X, we obtain 
the result 



1 



il/pof dk\a{k)\^k^'\ 



(4.22) 



V. EXACT SOLUTION 

The Fokker-Planck equation ( ^.Ij ) can be solved exactly for the Gaussian ensemble, 
which is the case of a parabolic potential V{E) = cE"^ (c is an arbitrary positive constant). 
The eigenfreauencies and eigenfunctions of the Fokker-Planck equation were constructed by 
Sutherland,lla by mapping it onto a Schrodinger equation. Here we use the same method 
to compute the correlation functions for /5 = 2, and compare with the asymptotic N oo 
results of Sec. |V[ 



A. Sutherland's method 



To map the Fokker-Planck equation ( |2.1| ) onto a Schrodinger equation we substitute 



P({K},r)=e-f^«^"»vp({Ej,r), 



(5.1) 



where W is given by Eq. ( p.2| ) with V{E) = cE"^. Sutherland^ used a different mapping 
(with P instead of P/2 in the exponent), but this one is more suitable for our purpose. 
Substitution of Eq. ( |0| ) into Eq. ( |2.1| ) yields for \Ef the equation 



9^ 
97 



1 A 



N 



f 


fdW^ 






\dEij 


I dE^ 



i=l ^-^i " / j,=l 

The expression between square brackets is evaluated as follows: 



(5.2) 
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= E E 7^^^ + E - 2ciV(iV - 1). (5.4) 

I \ ^ 3) t 

In the final equality we have used that for any three distinct indices i, j, k 

111111 , , 
I I = Q (5.5) 

Ei — Ej Ei — Ek Ej — Ei Ej — E^ E^ — Ei E^ — Ej 

so that the triple sum over k ^ i ^ j collapses to a double sum over i ^ j. Collecting 
results, we find that ^ satisfies a Schrodinger equation in imaginary time (r = it), 

- ^ = (^s - f/o)^, (5.6) 

«.^-^e|..^EE(^.^E-?, (-) 

U^ = N- + N{N -l)^. (5.8) 
7 ^7 

The Sutherland Hamiltonian Tis has an inverse-square interaction and a parabolic confin- 
ing potential. The interaction is attractive for (3 = 1 and repulsive for (3 = 4:. For (3 = 2 the 
interaction vanishes. Since exp(— /5W^) is a time-independent solution of the Fokker-Planck 
equation ( p.l| ), exp(— /5iy/2) is a time-independent solution of the Schrodinger equation ( p.6|) 
[in view of Eq. ([Ol) ]. Hence once has the eigenvalue equation 

Hse-^^ = f/oe-^'^. (5.9) 

For a particular ordering of the "coordinates" Ei, E2, . . . Ejy, the function \E'o oc exp{—pW/2) 
is an eigenfunction of the iV-fermion Hamiltonian Tis- Since it is nodeless, it is the ground 
state, at energy f/g. Anti-symmetrization yields the fermion ground-state wavefunctionlll'0 

M{En}) = Ce-|^«^"» n (5.10) 

with C a normalization constant. (Alternatively, we could work with the symmetric wave- 
function exp{—f3W/2), which is the ground state for hard-core bosons.) 

We obtain the A^-particle Green's function G{{En},T) of the Schrodinger equation ( ^.6|) 
from P{{En},T) by the similarity transformation 

G{{E„},t) = ef^«^"»P({E„},r)e-t^«^"». (5.11) 
For r > 0, the function G satisfies 
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8C 

- ^ = (^s - t/o)G, (5.12) 



in view of Eqs. ( |5.1| ) and ( |5.6| ). The initial condition is 



N 



G({E^},0) = l[S{E,-E% (5.13) 

i=l 

in view of Eq. ( ^^.6] ). Hence G is indeed a Green's function. In operator notation, 

G(r) = e-(«s-t/o)r_ ^5_^4^ 

We note that since the Fokker-Planck equation conserves the ordering of the levels 
Ei,E2,...En- for r > 0, we can write Eq. ( p.ll| ) equivalently in terms of the anti- 
symmetrized wavefunction ( |5.10 ), 

G({E„}, r) = %\{E^})P{{E^}, T)M{En})- (5.15) 

We are now ready to relate the equilibrium density-correlation function in Dyson's classi- 
cal Brownian-motion model to the ground-state density-correlation function in Sutherland's 
quantum many-body problem in imaginary time. In fact, we will see that the two correlation 
functions are identical. We define the ground-state expectation value {A)o of an operator A, 

/oo roc 
dE,--- dEj,%Am^. (5.16) 
-CXD J —OO 

The density operator is 



AT 



n 



{E)=Y.5{E-E,), (5.17a) 



n{E, t) = e^^^n{E)e-'^^\ (5. 17b) 

in the Schrodinger and Heisenberg picture, respectively. Combining Eqs. (|4.1|), (^^), ([4. 91), 
and (|5.14| )- (|5.17|) , one then finds 



(p(£;',r)p(E,0))eq = (n(E')e-(«^-^'^)^n(E)) 







(n(E',r)n(E,0))o. (5.18) 



Hence the density correlation function S{E, E' , X) defined in Eq. ( 4.10 ) is identical to 

S{E, E', X) = {n{E', X')n{E, 0))o - {n{E))o{n{E'))o. (5.19) 
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B. Gaussian unitary ensemble 

The significance of tlie formal relationship (|5.19| ) is that the quantum mechanical corre- 
lator on the right-hand-side can be computed exactly using the known excitation spectrum 
of the Sutherland Hamiltonian.lli The problem of computing the time-dependent correlation 
functions of Tio, was previously considered by Simons, Lee, and Altshuler, in connection with 
a microscopic theory of parametric correlations.0 We will return to their work in Sec. VII . 
The case (3 = 2 is particularly simple, since Tis is then the Hamiltonian of non-interacting 
fermions. This is the case of the Gaussian unitary ensemble. 

The single-particle eigenfunctions 4>p{E) and eigenvalues Ep of T^s are (cf. Eq. (|5.?1) with 
/? = 2): 

ME) = (2c/7r)i/^ i2Pp\)-'/\-''''%iEV2'c), (5.20) 
ep = {p+'^)u, p = 0,1,2,... (5.21) 
Lj = 2c/7. (5.22) 

The functions Hp(a;) are the Hermite polynomials. The density operator ( p. 17] ) becomes, in 
second quantization, 

oo 

n{E) = ME)ME)clc^, (5.23a) 

p,q=0 
oo 

n{E,T) = 0p(E)0,(E)e(^''-^«)-ctc^, (5.23b) 

p,q=0 

where cj^ and Cp are fermion creation and annihilation operators in state p. The average 
density in the A^-fermion ground state is 

N-l 

{niE))o^Pe^iE)=Y€iE). (5.24) 

p=0 

To compute the density fluctuations, we need the ground-state expectation value 

oo 

(n(E',r)n(E, 0))o = E ME')<PAE')ME)UE)e^''''~'''^^ ii'^A^qW (5-25) 
p',q',p,q=o 

The average of the product of four c's evaluates to 

(4'S'4s)o = ^q^W^iN - 1 - p')e{N - 1 - p) + 6p,,6p,,e{N - 1 - q)9{p - N), (5.26) 

where 5pq is the Kronecker delta and the function 9[x) equals 1 if x > and if a; < 0. 
Collecting results, we find for the density correlation function ( |5.19|) the formula 

oo N-l 

S{E, E', X) = Yi: ME)ME')UE)UE')e^'^-'''^''\ (5.27) 

p=N q=0 

The infinite series over p in Eq. (|5.27|) can be reduced to a finite sum by using an addition 
theorem for Hermite polynomials: 
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Go{E,E',t) = J2ME)<f)p{E')e-'^^ 

p=0 

c y/2 / c 

^TTsmhur) ^ Vsinhcur 



2EE' - {E' + E")coshuT 



(5.28) 



This is the famihar result for the (imaginary time) Green's function of a one-dimensional 
harmonic oscillator (with coordinate E, mass 7, and oscillator frequency u). Substitution 
into Eq. ( ^.271 ) yields 



N-l 

SiE,E',X) = GoiE,E',X') ^ 0,(E)0,(E')e^^''' 

q=Q 

N-l N-l 
p=0 q=0 



(5.29) 



For X = the function Gq becomes a delta function, so that Eq. (|5.29| ) reduces to 



N-l 



'N-l 



S{E, E\ 0) = 8{E - E') Y: <l>l{E) - E ME)ME') 



(5.30) 



p=0 



,p=o 



Eq. ( [5.30| ) was obtained by Mehtai using an approach known as the "method of orthogonal 
polymials". Our Eq. (|5.29|) extends this exact result to parametric correlations. 



C. Large- limit 

It is instructive to see how the result ( [4.14| ) of the asymptotic analysis in Sec. ^ follows 
(for (3 = 2) from the large-A^ limit of the exact result (|5.27|) for the Gaussian unitary 
ensemble. 

We wish to evaluate the density correlation function (|5.27|) in the limit N ^ 00, c — 0, 
while the product cN remains constant (to ensure a constant density of states, see below). 
Using the asymptotic form of the Hermite polynomials Hp for p ^ 1, one has for the 
eigenfunctions ( ^.201 ) the large-p expressions 

(l)2piE) = i-iy (2c/p7r2)V4e--^' cos(Ey8^), (5.31a) 
(p2p+i{E) = i-iy (2c/p7r2)i/^e-'=^' sin(E^). (5.31b) 

We need to compute the series 

N-l 5(^-1) 

E ME)ME'V''''' = E 02p(i?)02p(i?')e('^+^)"^' 

p=0 p=0 

|(7V-2) 

+ 02p+i(i?)02p+i(i?')e('^+^)"^', (5.32) 

p=0 
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in the limit — oo, c — ^ at constant cN . Note that c ^ imphes c<j ^ 0, in view of Eq. 
( ^.22| ). Combining Eqs. ( ^.31 ) and ( ^.32 ), and replacing the sum over p by an integral, we 
find in this limit 



p=0 

with the definitions 



n •'0 



po = -{cNY'\ 

TT 



a = Nu = 



27 ■ 



Similarly, 



p=N 



Po 



ds e 



COS {-Kp^iE - E')s). 



Substitution of Eq. into Eq. { ^2^ gives 

Peq(^) = Po, 



(5.33) 

(5.34) 
(5.35) 

(5.36) 
(5.37) 



justifying the identification ( 5.34 ). The limit N oo yields a uniform density of states 
in any fixed energy range. At finite A^, the density ppq(-E) vanishes for pol-^'l ^ 2N/ti, as 
follows from a more accurate evaluation of Eq. (|5.24|) .EI 

Substitution of Eqs. ( [5. 33] ) and ( |5.36| ) into Eqs. ( |5.27| ) gives an integral expression for 
the density correlation function S{E, E' , X) = S{E' — E,X), 

rl roc 

S{E,X)=pI ds ds' exp[aX'^{s'^ — s'^)]cos{TcpoEs) cos{TcpoEs'). (5.38) 
Jo Ji 

The Fourier transform S{k, X) = J dE S{E, X) exp{ikE) with respect to the energy incre- 
ment becomes 



S{k,X) 



Po 

e\k\ 



exp(-^^|/c|gmax) sinh(^^|/c|gn 



qmin = min 1, 



\k\ 



2vrpo 



= max 1, 



2vrpo 



(5.39) 
(5.40) 



The variable ^ was defined in Eq. ( 4.15|) . 

The result ( ^.391) holds in the large- A^ limit (at constant density of states) in any fixed k- 
range. We now further restrict ourselves to energy scales bigger than the mean level spacing 
A = pg \ i.e. to the range k <^ po. Eq. ( ^.39| ) then simplifies to 



S{k,X) 



\k\ 
2^ 



exp{-e\k\), 



(5.41) 



in agreement with Eq. ( [4.14| ) for f3 = 2. The correlation function S{k, X) is only appreciably 
different from zero if ^^|/c| ^ 1. Hence the restriction /c -C po on Eq. ( p.41| ) becomes irrelevant 
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if ^^po ^ 1- This implies that the asymptotic expressions for the correlation functions 
S{E,X) [and C{E,X)] of Sec. ^ hold for ^ oo in the energy range E > A for all X 
and in the parameter range X ^ for all E. If both E A and X ^ one can 

not use Eq. ( ^.41|) , but should use instead the full expressions ( ^.38|) or (|5.39|) . 



Once we have the density correlation function S{k,X), the current correlation function 
C{k,X) follows directly in view of the relation (|4.17|) . From Eq. ( |5.39|) one thus finds 



C{k,X) = ^^4|^°|3 exp(-^^|A:|gniax) (^^l^lgmin cosh(^^|fc|gmin)[3 + 4^^|A;|gmax] 

- sinh(f |fc|g^in)[3 + 3e\k\qm.. + 2ek\qL. + qLn)]) ■ (5-42) 



Eq. ( p.42| ) holds for N 00 and any k. For k <^ po it reduces to the asymptotic expression 
( CT ) of Sec. (with (3 = 2). 



VI. EXTENSION TO MULTIPLE PARAMETERS 

In this section we show how the Brownian-motion model of Sec. |I| can be extended 
to a parameter vector X = Xi,X2, . . -X^, relevant for a statistical description of the dis- 
persion relation of a d-dimensional crystalline lattice.0 The Brownian motion of the en- 
ergy levels En{X) then takes place in a fictitious world with multiple temporal dimensions 

f = Ti,T2, ...Td- 

We assume that any systematic drift in the energy levels is eliminated by a rescaling, so 
that 



d^E^iX) = 0. (6.1) 

(We abbreviate = d/dX^.) We also assume that the different parameters X^ are inde- 
pendent, that is to say 



dMX)duEn{X) = if /i 7^ I/. (6.2) 
Let r = coincide with X = 0. The initial condition on the distribution function P{{En}, t) 



IS 

N 



P{{E„},T = 0) = ll6{E^-E^), (6.3) 



i=l 



with Ef the eigenvalues of ?i(X = 0). For > the distribution function evolves according 
to the multiple-time-dimensional generalization of the Fokker-Planck equation ( |2.1| ), 

By comparing the initial average rate of change of the energy levels as a function of X and 



16 



we arrive as in Sec. |I| at the identifications 



(6.5) 
(6.6) 



(6.7) 
(6.8) 



2//?7m = (d.Ei)^. 

The Fokker-Planck equation (|6.4|) can now be reduced to a non-local diffusion equation 



as in Sec. IV 



(6.9) 



valid asymptotically for oo. For a constant density of states po the diffusion kernel 

becomes translationally invariant. Equation (|6.9| ) then has the fc-space solution 



Sp{k, t) = 6p{k, 0) exp 



-TTpo\k\ ^m/7a. 



(6.10) 



which implies for the density correlation function [cf. Eq. ( 4.14 )] 



^(fc,X) = ^exp 
up 



-7rpo\k\Y,Xl/% 



(6.11) 



Here we have also used the identification (|6.7] ). The i?-space correlation function becomes, 
upon inverse Fourier transformation, 



S{E,X) 



2n^/3 dE^ 
The current correlation function 



In ^^^[vrpolfclE^M/^.] 



(6.12) 



C^,iE,X,E\X') = Y.id^E,{X)){d,E,{X'))6iE - E,{X))5{E' - E,{X')) (6.13) 



is related to the density correlation function S{E, X, X') by 



dEdE' 



-C^,y{E,X,E',X') 



^2 



dX^dX, 



jS{E,X,E\X'). 



(6.14) 



Because of translational invariance, C^^{E, X , E' , X') = C^,j{E' — E,X' — X), 
S{E, X, E', X') = S{E' -E,X' -X). In A;-space Eq. ( |67[^ then takes the form 
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From Eqs. (|6.11| ) and (|6.15|) we find for tlie current correlation function the k and ii^-space 
expressions 



2po A 
1 92 



exp 



-npo\k\Y,Xl/-f^ 



ln{E' + [7rpo\k\Y^Xl/^, 



A=l 
l2 



(6.16) 
(6.17) 



For (i = 1 the correlation functions ( |6.12| ) and ( |6.17D reduce to the results ( |4.16| ) and 

{Km of Sec. m 



VII. COMPARISON WITH MICROSCOPIC THEORY 
A. Diagrammatic perturbation theory 

The asymptotic analysis of Sec. |V] yields the density and current correlation functions in 
the limit ^ cx) if E > A = 1/po for all X and if X > Xc = for all E. We will now 
show that the random-matrix theory (RMT) in this regime agrees with the diagrammatic 
perturbation theory of Szafer, Altshuler, and Simons.i'u 

When E 0, our result ( |4.19| ) for the current correlation function C{E,X) reduces to 

C(O..Y) = ^^ln|A-| 

^ itJf^O, (7.1) 



2' 



7r^f3X 

independent of the microscopic parameters po and 7. Eq. ( [7. 1| ) , obtained here from RMT, 
is precisely the universal correlator (|1.3|) which Szafer and Altshuleri derived from diagram- 
matic perturbation theory. 

At X = 0, the function C(0,X) according to Eq. (fTTI) has an integrable singularity 
consisting of a positive peak such that the integral over all X vanishes. This is a special 
case of the general sumrule 

/ dXC{E,X) = 0, (7.2) 
Jo 

which follows from Eq. ( [4.19| ) [cf. also Eq. ( |3.1CI|) ]. The peak of positive correlation has 
infinitesimal width in the limit E ^ 0. At non-zero E the peak has a finite width of order 
Xc(-E'/A)^/2, as illustrated in Fig. 2, where we have plotted C{E,X) from Eq. ( [4.19| ) for 
E = 0.1 A (dashed curve). 

As discussed in Sec. [V|, the asymptotic formula ( |7.1| ) becomes exact only for X 3> Xc. 
[Compare with the solid curve in Fig. 2, cornputed from the exact result ( ^.42| ).] Using the 
definition of the generalized Thouless energycl 
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Sc = ^-'El (7.3) 
and the relationship (p. 101) between 7 and E"^, one can write 



In Ref. the parameter X is the magnetic flux increment in units of h/e. Then is the 
conventional Thouless energy" Ec — hvpl/L'^, related to the conductance g (in units of 
e^/h) by (7 ~ Ec/A. The Aharonov-Bohm periodicity implies in this case the additional 
restriction X ^ 1 to Eq. ( |7.1| ) (which is compatible with the condition X ^ Xc because 
Xc — g~^^'^ -C 1 in the metallic regime). 

We have shown that the E ^ limit of C{E,X) obtained from RMT agrees with the 
microscopic theory. What about non-zero energy differences? This is most easily discussed 
in terms of the density correlation function S{E,X), to which C{E,X) is related via Eq. 
( CT )- Using ^ = X{^7ip£,y/^, we find that the result (CT) can be rewritten identically 
as 

SiE, X) = -^Re (iE + Itt(3£,X^) , (7.5) 

which agrees with the diagrammatic perturbation theoryll^ provided E <^ Ec- The deviation 
between RMT and the microscopic theory on energy scales greater than the Thouless energy 
Ec is well known from the work by Altshuler and Shklovskii on parameter-independent 
correlations.!! 

One can similarly show that the correlation functions for multiple parameters X^ {fi = 
1,2, ... d), obtained from RMT in Sec. |V|, agree with the results which Simons and Altshuleiil 
obtained by microscopic theory. In particular, we find from Eq. ( |6.17 ) that 



^ 2d - A 

E7.C,,(0,X) = -^[E^2/7M]"^ (7.6) 



in agreement with Ref. ^. The correlator ( |7.6|) is the multiple-parameter generalization of 
the universal correlator ( |1.3| ). Simons and Altshuler have discussed the physical origin of 
the different sign of the correlator for d < 2 and d > 2. 



B. Non-linear sigma model 

The restriction on the asymptotic analysis that either X ^ Xc or ^ A is removed 
by the exact solution of Sec. |V| for the Gaussian unitary ensemble (/5 = 2). The density 
correlation function in the limit X — 00 is given for this random-matrix ensemble by Eq. 
( p.38| ) in i?-space and by Eq. ( p.39| ) in fc-space. In terms of the generalized Thouless energy 
([7.3|), the -E-space expression can be written as 

S{E, X) = A-^ ds ds' exp 



2A 



cos{nsE/A)cos{TTs'E/A). (7.7) 
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This is precisely the result of the microscopic theory of Simons and AltshulerJ If either 
X S> Xc or A, Eq. ( |7.7|) reduces to Eq. ( |4.16| ) with [3 = 2. Simons and Altshuler were 
able to extend the microscopic theory to the regime X < Xc, -E ^ A, which is not obtainable 
by perturbation theory, by using a supersymmetry formulation followed by a mapping onto 
a non-linear sigma model.H As emphasized by these authors, it is quite remarkable that the 
microscopic parameters enter only via the quantities A and £ci so that a rescaling of the E 
and X variables maps all density correlation functions onto a single universal function. 

Simons and Altshuler have also computed the small E and X behavior of S{E,X) from 
the microscopic theory in the presence of time-reversal symmetry, i.e. for (3 = 1. Again, they 
used a mapping onto a non-linear sigma model to go beyond perturbation theory. We have 
no RMT result for the small E and X behavior in this case, which would correspond to the 
orthogonal ensemble. The case /3 = 4 of strong spin-orbit scattering (symplectic ensemble 
in RMT) has not yet been treated by microscopic theory, and only in the asymptotic limit 
by RMT. 

In Ref. |2y, Simons, Lee, and Altshuler have argued convincingly (although not com- 
pletely proven) that the density correlation functions of the non-linear sigma model and the 
Sutherland Hamiltonian are equivalent for (3 = 1,2, and 4. In the present paper, in Sec. 0, 
we have proven the equivalence of the density correlation functions of Dyson's Brownian- 
motion model and the Sutherland Hamiltonian. Taken together, this is evidence for the 
complete equivalence of the non-linear sigma model and RMT, although the cases (3 = 1 
and (3 = 4 still lack a complete proof. 



C. Conclusion 

We have studied the response to an external perturbation of the energy levels of a quan- 
tum mechanical system, by means of the Brownian-motion model introduced by Dyson 
in the theory of random matrices. Our results for the energy and parameter-dependent 
level-density and current-density correlation functions S{E,X) and C{E,X) agree with the 
microscopic theory for a disordered metallic particle, for energy scales below the Thouless 
energy E^.. This establishes the validity of Dyson's basic assumption, that parametric cor- 
relations are dominated by level repulsion and therefore solely dependent on the symmetry 
of the Hamiltonian. 

It is likely that the approach developed in this paper can also be used to describe para- 
metric correlations in random transmission matrices. The analogue of level repulsion for 
the transmission eigenvalues is known,!! and leads to a pair correlation function K(T, T') 
which differs from Eq. ( p.. 61 ) for K{E, E') but has the same universal /3-dependence.li3 This 



suggests that the analogue of the universal correlator ( |1.3|) exists as well for the transmission 



eigenvalues, with obvious implications for the conductance of a mesoscopic system. 
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APPENDIX: DERIVATION OF EQUATION 



For completeness, we present here Dyson's derivationEJ of the non-hnear diffusion equa- 
tion ( [4 .21 ) from the Fokker-Planck equation ( |2.1D , in the hmit oo. 

We multiply Eq. by 6{E — Ei), integrate over Ei, E2, . . . Ejy, and sum over i. The 
result is 



(9 /'°° 
(3-^ p{E,t)+ dE,--- dENP{{E^],r)Y,6{E-E, 

Oil/ J —00 J —00 



dW 



dEi 



(Al) 



where p{E,t) is defined in Eq. (|4.1|). Substitution of the definition ( |2.2| ) of W{{En}) into 
Eq. (^) leads to 



P-'g^piE, r) + p(E, r)—V{E) - V j jE 



>P 2{E,E',t) 
E-E' 



(A2) 



where V J indicates the principal value of the integral. The pair density p2iE, E', r) is defined 

by 

/oo roo 
dE,--- dENP{{Er,},T)J26{E-E,)6{E'-E^). (A3) 
-00 J —00 ■/ ■ 

The pair density is symmetric in the energy arguments, P2{E,E',t) = P2{E',E,t), and 
satisfies the normalization 



/oo 
dE'p2{E,E',T) = {N-l)p{E,T). 
-00 



(A4) 



Following Ref. |T3| we decompose the pair density into a correlated and an uncorrelated 
part, 

P2{E, E', t) = p{E, T)p{E', t) [1 - y{E, E', r)]. (A5) 
The function y{E, E', r) = y{E', E, r) is symmetric in E and E', and satisfies 

/oo 
dE'y{E,E',T)p{E',T) = l, (A6) 
-00 

in view of the normalization ([A4|) . Substitution of the definition (|A5|) into Eq. (|A2|) leads 

to 



^"'dE^^^' ^) + PiE, ^)Q^\yiE) + U{E, r)] + p{E, t)I{E, r) 



, (A7) 



21 



with the definitions 

/(E,r) =p/_^rfEV(E^r)^^|^, (A8) 

/oo 
dE' p{E',r) In \E - E'\. (A9) 
-oo 

Eq. ( 1^) is still exact. To introduce the approximation we need one further piece of 
notation. We re-express the function y{E, E', r) in terms of the sum and difference variables 
t = ^iE + E') and s = E'-E: 

y{E, E', r) = Y{l{E + E'), E' -E,t) = Y{t, s, r). (AlO) 

The function Y{t, s, r) = Y{t, —s, r) is even in s. The normalization ( [A 61 ) becomes 

/oo 
dsY{E + ls,s,T)p{E + s,T) = l. (All) 
-oo 

Similarly, the integral (|A8|) takes the form 

/oo 
dsY{E + ls,s,T)p{E + s,T)s-\ (A12) 
-oo 

By substituting the Taylor expansions 

Y{E + Is, s, t) = Y{E, s, t) + \s^Y{E, s, r) + ■ • • , (A13) 
p{E + 5, r) = p{E, r) + s^p{E, r) + ■ ■ • , (A14) 

into Eq. ( |A12| ), we obtain an expansion of I{E,t) in higher and higher moments Yp{E,T) 
of Y{E, s, r) with respect to s, 

/oo 
dsY{E,s,T)sP. (A15) 
-oo 

Because of the symmetry Y{t, s, r) = Y{t, —s, r) only even moments contribute [Yp{E, r) = 
for p odd]. Following DysonjB we neglect the second and higher moments. An order of 
magnitude estimate suggests that the error involved in neglecting Yp for p > 2 is of order 
N~^. Dyson argues that the error is actually of order A^~^lnA^, by comparison with exact 
results for the distribution of the spacing of eigenvalues. 

Since Y_i and Yi are identically zero, only Yq contributes to I{E, r) to second order. 
Substitution of the Taylor expansions (|A13| ) and (|A14|) into Eq. ( |A12| ) yields 



I{E, r) = -lp(E, r)^Yo{E, r) - Y,{E, r)—p{E, r). (A16) 



Similarly, substitution of the Taylor expansions into Eq. ( [A 111) yields 

p{E,t)Y,{E,t) = 1. (A17) 
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Combining Eqs. ([A 161) and ( |A17|) we find 



(A18) 



Hence Eq. 
d 



This is Eq. 



takes the form 
d 



d_ 
'dE 



dE 



V{E) - / dE'p{E\ r) \yi\E-E'\ + 



2-£ 

2/3 



lnp{E,T) 



(A19) 



except for the final term, proportional to (2 — f3)/2f3. As noted in Ref. 
T3| , this term is of order In and can be neglected relative to the other terms, which are of 
order A^. Dropping that term, we obtain Eq. ( [4.2|) . 
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